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The paper, authored by J. A. S. Lima et al, was published in Phys. Rev. E in 2000 has discussed 
the dispersion relation and Landau damping of Langmuir wave in the context of the nonextensive 
statistics proposed by Tsallis. It has been cited by many authors because the dispersion relation in 
Tsallis formalism present a good fit to the experimental data when q < 1, while the classical result 
based on Maxwellian distribution only provides a crude description. However, the results obtained 
in this paper are problematic. In this comments on the paper we shall derive the correct analytic 
formulas both for the dispersion relation and Landau damping in Tsallis formalism. We hope that 
this comments will be useful in providing the correct results. 
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I. INTRODUCTION 

Over the last few years, it has been proven that sys- 
tems which present long-range interactions, long-time 
memory, fractality of the corresponding space-time, or 
intrinsic inhomogeneity are intractable within the con- 
ventional Boltzmann-Gibbs statistics [1, 2]. So there has 
been an increasing focus on the new statistical approach, 
i.e. nonextensive statistical mechanics (NSM), in recent 
years. For q ^ 1, it gives power-law distribution and only 
when the parameter q — > 1 Maxwellian distribution is re- 
covered [3]. NSM has been successfully applied to stellar 
polytropes [4], two dimensional Euler and drift turbu- 
lence in a pure electron plasma column [5], as well as to 
the peculiar velocity function of galaxies clusters [6] . In 
particular, Liu et al. [7] showed a reasonable indication 
for the non-Maxwellian velocity distribution from plasma 
experiments. 

Dispersion relations are fundamental and important 
for studying the wave in the plasma. According to the 
dispersion relations, we can study the problem of instabil- 
ity, propagation, refraction and absorption of the plasma 
wave. Recently, there has been a great deal of inter- 
est in studying the dispersion property in plasmas in the 
context of the nonextensive statistics. The paper au- 
thored by J. A. S. Lima et al [8], has a considerable hold 
over the field of plasma physics, have studied the disper- 
sion relation of Langmuir wave based on q-distribution, 
the results show that nonextensive formalism presents 
a good fit to the experimental data, while the stan- 
dard Maxwellian distribution only provides a crude de- 
scription. However, the results obtained in this paper 
are problematic. The reason is that the author inves- 
tigates the propagation of electrostatic waves by using 



the one-dimensional equilibrium distribution. But the 
equilibrium distribution in the dielectric function should 
be marginal distribution (quasi one-dimensional distribu- 
tion), which is different with one-dimensional distribu- 
tion in the context of nonextensive statistics unlike the 
classical Boltzmann-Gibbs statistics. In this comments 
on the paper we will derive the correct analytic formulas 
both for the dispersion relation and Landau damping in 
detail with Tsallis formalism. It is our hope that the dis- 
cussion here will be useful in the field of plasma physics. 

The paper is organized as follows. In Section II we 
briefly introduce the nonextensive distribution function. 
The generalized dispersion relation and Landau damping 
for Langmuir wave are obtained in Section III. Finally, 
the summary is given in Section IV. 



II. NONEXTENSIVE DISTRIBUTION 
FUNCTION 

First let us recall some basic facts about Tsallis statis- 
tics. In Tsallis statistics, the entropy has the form [3] 
of 
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where kg is the Boltzmann constant, q is a parameter 
quantifying the degree of nonextensivity, pi is the proba- 
bility of the ith microstate. The B-G entropy is recovered 
in the limit q — > 1. The basic property of Tsallis entropy 
is the nonadditivity or nonextensivity for q ^ I. For ex- 
ample, for two systems A and B, the rule of composition 
[3] reads 

S q (A + B)= S q (A) + S q (B) + (1 - q)S q (A)S q (B). (2) 

In the nonextensive description, the three-dimensional 
equilibrium distribution function can be written as [9] 
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f q (p)=A q [l-(q-l) 



according to the normalizing condition 
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where Eq. (fTU|) has been calculated using the integral for- 
mula [10, p325], that is, 
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with |arg/3| < 7r,p > 0,0 < i?e/i < pRev, here require 
q > 3/5. -B is the Beta function, the relation of Beta 
function and Gamma function is [10, p909] 
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where p, «t = \JksT Jm, fee, T, to, and no denotes, 
respectively, the momentum of particles, thermal speed, 
Boltzmann constant, temperature of particles, mass of 
particles and particle number density. As one may check, 
for q < 1/3, the q-distribution is unnormalizable. For 
1/3 < q < 1, the momentum of the particles can take 
any value. For q > 1, the distribution function Eq.Q 
exhibits a cutoff on the maximum value allowed for the 
momentum of the particles, which is given by 
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For q > 1, 
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We see that in the limit q — > 1, p max goes to infinity and 
Eq.Q reduces to the Boltzmann distribution function 
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In order to define the temperature of the system which 
is described by the nonextensive distribution, we will cal- 
culate the average kinetic energy below. For 1/3 < q < 1, 
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where Eq. fTB")) has been calculated using the transfor- 
mation t = (q — 1) p 2 j2rv?v\ and integral formula [10, 
p324], that is, 



3 



c"- 1 fl 



dx = 



>6H 



(14) 



with i?e/x > Q,Rev > 0, A > 0. So the average kinetic 
energy can be expressed as 
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then the integral in Eq. (fl"9]) over p y is 



°° 2m 2 v 2 T + (l-q)(p 2 x+ pl) (l- g ) , 
„ 1 2m2«2 + 2m 2 4 PyJ ^ 



where T 9 = 2T/(5q — 3) is the physical temperature of 
the nonextensive system. We see that in the limit q — > 
1, T q= i — T and the average kinetic energy reduces to 
(E q= i) = 3nokBT/2, which is the standard result in B-G 
statistics. 



III. THE GENERALIZED DISPERSION 
RELATION AND LANDAU DAMPING 

For the longitudinal wave propagating in an unmag- 
netized, collisionless, isotropic plasma, the longitudinal 
dielectric function of electron can be written as [11] 
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we consider the direction of wave vector k to be along 
x-axis, Eq. p^)) becomes 



{2m 2 v 2 T + (l-q)(pl+pl)y^ + \ (20) 

where Eq.(|20|) has been calculated using the integral for- 
mula (fTTj) . and substituting Eq. ([2U)) into Eq. (ITO|) . accord- 
ing to the same method we can calculate the integral over 
p z . Finally Eq. (fl"9l) becomes 
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where e is the electron charge, id comes from Landau 
rules(<5 — > + ) [12]. Note that f q (p x ) is the marginal dis- 
tribution in the nonextensive framework, which is given 
by 
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then the integral in Eq.(|22|) over becomes 
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Next we will derive the expression of the marginal dis- 
tribution. Substituting Eq.© into Eq.(|T8]). for 3/5 < 
Q < 1, we obtain 
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unlike the classical Boltzmann-Gibbs statistics. It is the 
reason why the results obtained by J. A. S. Lima et al, 
are problematic. 

Substituting the marginal distribution Eqs. (|21l) and 
(|24|) into the dielectric function Eq. (fT7|) , we obtain 
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where w pe = ■J Airnoe 2 / m is the plasma frequency, x 

is the dimensionless parameter, namely, x = Loj\[2kvT- 
Z q ix) is the generalized plasma dispersion function in 
the context of Tsallis statistics, 
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where Eo. ([23[) has been calculated using the transfor- 
mation t = (q - \) p 2 J {2m 2 v\ - (q - 1) (p 2 + p 2 ) ) and 
integral formula (jT4[) . then substituting Eg. (123)) into 
Eq. (|22|) , according to the same method we can calculate 
the integral over p z . Finally Eg. (f2"2"]l becomes 
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(29) 

where £ = v x /v2vt, in the limit g — > 1, it is reduced to 
the standard form in B-G statistics [11] 
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where p denotes the principal value, then the generalized 
plasma dispersion function Eq. (|29j) can be written as 
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Obviously, the marginal distribution Eqs. (|21|) and (|24|) 
are different with one-dimensional distribution in the 
context of nonextensive statistics [9] 
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when lo kvx, namely x 3> 1, the real part of Eq. ([32 
becomes 
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When bj <C kvx, namely introducing the trans- 

formation £ = r/ + x, then the real part of Eq. ([3"2"|) can be 
written as 



q < 0.85, obviously, it should be 0.82 < q < 0.91 based 
on the correct result. 

Next we will derive the expression of Landau damping. 
The Landau damping rate can be written as [11] 
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we obtain the generalized Landau damping as 
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combined with Eq. (f3T)| , 



7fe = -L q J -UJ pe — 



which can be used in investigating the low-frequency 
wave, such as the ion acoustic waves. It should be noted 
that the process is not pinpoint in Eq. ([35"T) . the real part 
should be a very small quantity, which may be obtained 
by numerical method. However, the small quantity can 
be neglect when Eq. ([32| Substituted into the dielectric 
function Eq.(|28]). 

Substituting Eq.(|34| into the dielectric function 
Eq.([2"8]), according to the longitudinal dispersion rela- 
tion Ree l k = 0, thus the generalized dispersion relation of 
Langmuir wave is obtained, 
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where kd = uj pe /vT is the electronic Debye wave number. 
In the limit q —> 1, Eq. (|40)) reduces to 
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(41) 



which is the classical Landau expression for the damping 
decrement in the framework of B-G statistics [11]. 



IV. SUMMARY 



pe 
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where vr q — \/kBT q /m is the physical thermal speed, 
T q is the physical temperature defined in Section 2. As 
expected, in the limit q — > 1, Eq.(|3"7]) reduces to 



uj 2 + 3k 2 v T , 



(38) 



being the standard result in B-G statistics [11]. Thus the 
dispersion relation of Langmuir wave obtained by J. A. 
S. Lima et al [81 



3q - 1 



-3fcV 



is problematic. According to Fig.l in their paper, we 
can see that the dispersion relation for Tsallis formalism 
presents a good fit to the experimental data when 0.7 < 



In this comments, we have discussed the dispersion 
property and Landau damping of Langmuir wave in an 
unmagnetized, collisionless, isotropic plasma with the 
nonextensive distribution in Tsallis statistics. The cor- 
rect generalized dispersion relation and Landau damping 
are obtained. In the limiting case (q —> 1) the classical 
results based on the B-G statistics are recovered. It is 
our hope that the discussion here will serve as a useful 
introduction to the field of plasma physics. 
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